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Preface 


This  thesis  is  the  result  of  my  efforts  to  analyze  the 
structure  of  several  amorphous  alloys  through  the  use  of 
Mosshauer  spectroscopy.  There  were  four  objectives  in  this 
study;  (1)  to  obtain  the  Mossbauer  spectra  of  each  sample  at 
room  and  at  liquid  nitrogen  temperatures;  (2)  to  analyze  the 
obtained  spectra,  relating  the  results  to  the  structure  of 
the  materials;  (3)  to  modify  the  existing  Mossbauer  equipment 
to  simplify  the  use  and  adjustment;  and  (^)  to  simplify  the 
computer  program  used  in  data  analysis.  This  thesis  will 
describe  the  theory  of  the  amorphous  structure,  and  the  model 
used  to  describe  it,  the  equipment  and  experimental  methods 
used,  the  analysis  of  the  data,  the  results,  and  the  con- 
clusions and  recommendations. 

I would  like  to  acknowledge,  and  thank,  several  people 
who  were  most  helpful  in  this  project.  Dr.  Harold  Gegel,  of 
the  Air  Force  Materials  Laboratory,  for  supplying  the  materials 
used  in  this  study;  Dr.  Jon  Spijkerman,  of  Ranger  Electronics, 
for  his  material  and  personal  support;  and,  finally.  Dr. 

George  John,  without  whose  help  and  guidance  this  project 
could  not  have  been  completed. 

I also  wish  to  express  my  thanks  to  my  wife,  Denise,  for 
her  unwaiver ing  support  and  patience  in  this  undertaking. 
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Abstract 

Several  models  have  been  used  to  describe  the  structure 
of  iron-phosphorous-carbon  amorphous  alloys.  In  this  study 
a twelve  site  substitutional  solid  model  is  used  to  describe 
the  Mossbauer  spectra  obtained  from  six  different  metallic 
amorphous  alloys.  Spectra  were  obtained  at  room  and  at 
liquid  nitrogen  temperatures  for  the  parent  material,  as  well 
as  for  alloys  containing  various  percentages  of  boron  or 
silicon.  Each  spectrum  is  fit  with  a sum  of  Lorenzian  pro- 
files describing  the  twelve  site  model.  The  results  indicate 
that  the  model  is  inadequate  in  fully  describing  the  structure 
of  the  various  materials.  The  results  also  indicate  a dif- 
ference in  structure  between  the  materials,  and,  in  two 
cases,  significant  changes  in  structure  with  decreased 
temperature . 
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AN  ANALYSIS  OF  METALLIC  GLASSES 
BY  MOSSBAUER  SPECTROSCOPY 

I . Introduction 

Mossbauer  spectroscopy  is  used  as  an  analytical  tool  in 
disciplines  as  far  ranging  as  biology  and  lunar  seology. 

The  wide  use  of  the  technique  stems  from  the  method's  abil- 
ity to  provide  information  on  the  environs  of  an  individual 
nucleus.  The  analysis  of  a Mossbauer  spectrum  can  give 
information  on  the  chemical,  electrical,  and  magnetic  envi- 
ronment of  the  Mossbauer-isotope  nucleus.  Because  of  this 
unique  ability,  Mossbauer  spectroscopy  has  become  a major 
tool  in  physical  metallurgy.  The  analyst  hopes  to  extrap- 
olate from  the  various  Mossbauer  parameters  obtained,  infor- 
;•  mation  descriptive  of  the  macroscopic  material. 

Background 

Amorphous  metallic  alloys,  sometimes  called  the  me- 
tallic glasses , exhibit  several  properties  which  recommend 
their  use  as  structural  materials.  Among  these  are  flex- 
ibility and  high  strength  (as  high  as  2 to  4 x 10-^  psi)  — 
both  highly  desirable  characteristics  (Ref  1;563)-  It  is 
these  properties  that  have  attracted  the  attention  of  the 
Air  Force  Materials  Laboratory  to  the  amorphous  metallic 
alloys  for  possible  use  in  Air  Force  weapons  systems. 

The  Materials  Laboratory  has  begun  a study  of  these 
alloys,  and  as  a part  of  that  study,  is  examining  a series 
of  amorphous  compounds  manufactured  by  the  Battelle 
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Laboratories  at  Columbus,  Ohio.  These  materials  consist  of  a 
parent  alloy  of  iron-phosphorous-carbon  and  a series  of 
daughter  alloys  consisting  of  the  parent  alloy  with  added,  or 
substitutional,  boron  or  silicon.  The  presence  of  iron  in 
each  of  the  samples  permits  the  use  of  Mossbauer  spectro- 
scopy with  a cobalt-57  source  without  modification  to  the 
sample . 

To  fully  predict  the  response  of  these  materials  to  a 
given  environment,  a knowledge  of  their  internal  structure  is 
desirable.  However,  the  structure  of  the  metallic  glassy  al- 
loys has  not  been  fully  described.  Various  models  describing 
the  structure  have  been  proposed  (for  example.  Ref  2,3);  how- 
ever, none  of  these  models  has  been  totally  successful,  al- 
though each  of  the  models  does  provide  a means  to  describe  a 
property  or  reaction  of  an  amorphous  alloy  to  a given  stimu- 
lus. This  study  uses  one  of  the  many  models,  primarily  for 
its  simplicity,  and  does  not  address  its  validity  except  in 
terms  of  the  results  of  this  study. 

Purpose 

The  primary  purpose  of  this  study  was  to  determine  what 
can  be  learned  of  the  structure  of  amorphous  materials 
through  the  use  of  Mossbauer  spectroscopy.  Within  this  broad 
purpose,  the  specific  objectives  were  (1)  to  obtain  the 
Mossbauer  spectra  of  the  various  compounds  at  room  and  at 
liquid  nitrogen  temperatures,  and  (2)  to  begin  an  analysis 
relating  the  data  obtained  to  the  structure  of  the  compound. 

Secondary  aims  were  to  modify  the  existing  Mossbauer 
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equipment  in  order  to  improve  spectrum  resolution  and  ease 
of  adjustment,  and  to  simplify  and  upgrade  the  computer  pro- 
gram used  in  data  analysis  and  curve  fitting. 

Overview 

The  theory  and  use  of  the  Mosshauer  effect  using  a co-  t 

balt-57  source  and  an  absorber  containing  iron-57  are  well 
documented  (e.g.  Ref  4,5)  and  will  not  be  discussed  in  this 
study.  The  completion  of  this  study  required  the  selection  i 


of  a model  of  the  internal  structure  of  the  amorphous  mater- 
ials. That  model  and  the  associated  theory  are  discussed  in 
Chapter  IT.  The  Mossbauer  equipment  used  and  the  modifica- 
tions made  to  the  equipment,  as  well  as  the  method  of  data 
collection,  are  discussed  in  Chapter  III.  Chapter  IV  covers 
data  processing  and  analysis,  along  with  the  significant 
problems  encountered  in  describing  the  experimentally-obtain- 
ed spectra.  The  remaining  chapters  discuss  results,  conclu- 
sions, and  recommendations. 
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II.  Theory  and  F/lodel 

Various  models  have  been  used  to  describe  the  structure 
of  amorphous  materials.  This  descussion,  however,  will  be 
limited  to  several  which  are  applicable  to  the  glassy  metal- 
lic alloys.  In  particular,  this  study  is  concerned  with 
those  models  which  can  be  used  to  describe  iron-phosphorous- 
carbon  amorphous  materials. 

In  an  amorphous  material,  the  atoms  have  coagulated  into 
an  extraordinarily  disordered  state  rather  than  the  normal 
crystalline  state.  This  disordered  state  is  normally  achiev- 
ed by  the  rapid  cooling  of  the  molten  state  by  splat  quench- 
ing, spin  cooling,  or  vapor  deposition  (Ref  6:231).  Polk 
and  others  have  described  the  resulting  material  as  a solid- 
ified supercooled  liquid  (Ref  7:k86) . The  diffuse  halo  dif- 
fraction patterns  obtained  by  Fujita  and  others  from  various 
amorphous  alloys  exhibit  random  distribution  patterns  clearly 
different  from  that  of  the  liquid  from  which  the  alloys  were 
made  (Ref  6:231).  In  both  cases  their  results  implied  a ran- 
dom distribution  function. 

Polk  proposed  a structural  model  consisting  of  "a  ran- 
dom packing  of  the  metal  atoms,  which  would  have  an  arrange- 
ment similar  to  the  hard  sphere  packing  studied  by  Bernal, 
with  the  metalloid  occupying  the  larger  holes  inherent  in 
such  a structure."  The  hard  sphere  packing  referred  to  by 
Polk,  and  described  by  Bernal,  has  idealized  holes  in  the 
form  of  tetrahedrons,  octahedrons , capped  trigonal  prisms, 
and  several  even  more  complex  geometries  (Ref  7:^8?), 
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Fujita,  in  studying  an  amorphous  alloy  similar  to  that  of 
this  report,  modeled  the  structure  as  a body  centered  cubic 
with  the  metalloid  atoms  substitutional  in  the  cubic  struc- 
ture, as  well  as  existing  interstitially  (Ref  6:232).  Tsuei, 
in  an  earlier  work  done  in  1968,  also  assumed  the  body  center- 
ed cubic,  but  treated  only  the  eight  nearest  neighbors,  as 
opposed  to  Fujita's  treatment  of  a total  of  nearest  and 
second  nearest  neighbors  (Ref  8:605). 

Model  Selection 

A basic  question,  then,  in  the  selection  of  a model  to 
use  in  the  analysis  of  the  experimental  data,  is  the  number 
of  first  and  second  neighbor  atoms  to  be  considered.  The 
Mos.sbauer  spectrum  obtained  from  a given  material  is  a result 
of  the  total  interaction  of  all  of  these  neighbors  with  the 
absorber  atom;  and  the  model  must  account  for  each  individual 
effect . 

Tsuei's  study  of  FegQP]^2 .5^7 .5  (where  the  subscripts 
indicate  atomic  percent),  as  noted  before,  assumed  eight 
nearest  neighbors,  but  used  only  the  five  most  prominent  com- 
binations of  iron  and  non- iron  atoms  in  his  computer  curve 
fitting.  Here,  and  throughout  this  study,  the  most  prominent 
combinations  are  those  which  have  the  greatest  effect  on  the 
Mossbauer  absorber  atom  and  will  be  referred  to  as  the  most 
prominent  iron  "sites".  The  results  obtained  by  Tsuei  approx- 
imated a binomial  distribution  of  iron  and  non- iron  atoms  in 
the  nearest  neighbor  shell  (Ref  8)  . 

Complicating  Tsuei's  results  was  the  determination  by 
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Lin  that  a similar  alloy  could  be  described  by  13*6  nearest 
neighbors  consisting  of  a random  distribution  of  iron  and  non- 
iron atoms.  Fujita,  in  discussing  these  results,  reports 
that  he  was  able  to  analyze  the  Mossbauer  patterns  for  iron- 
phosphorous-carbon  alloys  by  assuming  a "nearly  perfect  body 
centered  cubic  consisting  of  very  small  crystallites."  This 
assumption  lead  to  the  consideration  of  eight  first  neighbor 
and  six  second  neighbor  sites  (Ref  6:232).  While  Fujita's 
results  work  well  in  fitting  the  experimental  data,  two  prob- 
lems exist  in  the  application  of  hio  scheme. 

The  first  problem  is  expense.  The  equation  used  to  de- 
scribe the  spectrum  becomes 
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l+Q(X-X-6^.n  -6  m )/r']‘ 


(1) 


The  meaning  of  these  terms  will  be  discussed  later.  What  is 
significant  here  is  that  an  iterative  best-fit  computer  scheme 
would  require  in  excess  of  1.5  million  operations  per  itera- 
tion of  the  22  possible  variables.  Unless  the  values  of  the 
unknowns  are  known  to  a reasonable  degree  prior  to  iteration, 
this  scheme  would  require  a great  deal  of  computer  time. 

The  second  problem  is  interpretation  of  the  results. 

The  calculated  spectrum  would  represent  the  sum  of  288 
Lorenzian  peaks.  Intuitively,  one  feels  any  spectrum  could 
be  fit  using  this  many  peaks. 


Simplified  Model . Fujita  also  describes  a simpler  model 
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applicable  to  substitutional  solid  solutions.  In  this  model, 
the  effects  of  the  second  neighbors  are  assumed  to  be  neg- 
ligible; and  a random  distribution  of  solute  atoms  are  assum- 
ed to  exist  in  a closely  packed  substitutional  solution.  The 
atoms  are  treated  as  hard  spheres  of  equal  diameter  packed  to 
near  maximum  density.  The  resulting  structure  yields  12  near- 
est neighbor  atoms.  Fujita  develops  the  probability  of  find- 
ing a given  combination  of  atoms , assuming  that  only  two  types 
are  present.  The  probability  is 

P(nA,inB)=j^2Cn^^(^”°)"^  » n+m=12  (2) 

where  n is  the  number  of  atoms  of  species  A,  m is  the  number 
of  species  B,  C is  the  binomial  coefficient  for  12  items  taken 
n at  a time,  and  c is  the  concentration  of  species  A.  These 
definitions  apply  to  Eq  (1)  as  well.  With  the  probabilities 
defined  in  this  manner,  the  intensity  at  the  calculated  spec- 
trum can  be  described  by 

" A „ n.,  .12-n  i , . 

I(x)-  2,T“12^n°  (I-"")  l+[(x-Xo-6^)/r]^  (3) 

n=0 


where  A is  the  total  absorption  intensity,  Xq  is  the  peak 
position  with  no  foreign  neighbors,  represents  the  line 
shift  due  to  one  foreign  neighbor,  and  r is  the  half  width 
at  half  maximum  of  the  calculated  Lorenzian.  Using  Eq  {3)» 
the  calculated  spectrum  is  the  sum  of  12  Lorenzian  peaks. 

The  contribution  from  each  site  is  a single  peak  spec- 
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trum,  which  is  not  what  is  expected  from  an  iron-rich 
alloy,  nor  is  it  consistent  with  the  results  of  previous 
work.  Several  final  assumptions  and  modifications  are  re- 
quired . 

Final  Model . The  first  modification  to  Eq  (3)  is  based 
on  the  assumption  that  the  contribution  from  each  site  will 
be  a six  peak  spectrum  similar  in  form  to  that  produced  by 
natural  iron,  and  describable  by  the  simplified  Hamiltonian, 
as  developed  by  Wertheim  (Ref  9572-83).  This  assumption  is 
based  on  the  results  of  Tsuei  and  Fujita,  and  is  considered 
to  be  a valid  one.  The  remaining  assumptions  determine  which 
terms  of  the  final  expression  []Eq  (4)]  shall  be  considered 
variables . 

The  first  assumption  is  that  each  of  the  sites  shall 
have  its  peak  intensities  in  the  same  ratio,  but  that  ratio 
will  not  be  fixed.  If  the  spins  in  the  material  were  ran- 
domly oriented,  a ratio  of  352.‘1: 1:2:3  for  the  six  peaks  con- 
tributed by  each  site  would  be  expected.  Tsuei 's  results 
indicated  a peak  ratio  other  than  this,  and  therefore  the 
simplified  model  of  this  study  was  designed  to  allow  arbit- 
rary variation  of  the  peak  ratios.  The  arbitrary  decision 
to  require  the  peak  heights  of  each  site  to  be  in  the  same 
ratio  was  made  to  simplify  the  model  by  reducing  the  number 
of  variables  for  peak  height  from  72  to  6.  Note  that  this 
scheme  does  not  require  symmetrical  peaks  to  be  of  the  same 
intensity. 

The  second  assumption  is  that  linewidth  is  an  arbitrary 
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variable.  The  model  was  designed  to  allow  three  values  of 
linewidth  to  vary  without  restriction.  Each  of  the  three 
represents  a pair  of  symmetrical  peaks  (e.g.  peak  1 and  peak 
6)  for  a given  combination.  Each  of  the  combinations  is  re- 
quired to  have  the  same  linewidths , again  an  arbitrary  de- 
cision to  increase  simplicity  and  decrease  the  number  of 
variables.  While  the  assumption  does  simplify  the  model,  it 
is  unreasonable  in  the  sense  that  the  model  allows  two  sym- 
metrical peaks  to  vary  independently  in  amplitude,  yet  re- 
quires each  to  have  the  same  full  width  at  half  maximum.  If 
the  experimental  data  is  in  the  form  of  an  asymmetrical  pat- 
tern, the  values  found  by  the  least-squares  minimization  pro- 
gram will  represent  only  a weighted  average,  and  not  the  best 
possible  fit. 

The  third  assumption  is  that  the  contributions  from  the 
six  least  prominent  combinations  could  be  neglected.  This 
assumption  was  made  to  reduce  execution  time  of  the  curve- 
fitting program,  and  should  introduce  minimal  error  into  the 
results.  The  weighting  function  for  the  sixth  most  prominent 
combination  is  0.0532,  and  the  sum  of  the  contributions  from 
the  six  remaining  combinations  is  O.OI93. 

The  final  assumption  is  that  the  constituent  atoms  of  the 
materials  to  be  analyzed  can  be  treated  as  iron  and  non- iron. 
This  implies  that  the  shift  in  energy  levels  of  the  absorber 
nucleus  is  independent  of  the  type  of  atom  interacting  with 
the  absorber  — which  is  contradictory  to  previous  results 
with  Mossbauer  spectroscopy,  but  does  allow  use  of  the  bi- 
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nomial  distribution.  The  validity  of  this  assumption  is 
further  discussed  in  the  Conclusions. 

The  final  form  of  the  model  is  then 


12  6 


n^m.i  ^ 1+  [ (x-Xq- 6Q,n)/  r] 


12-n 


(^) 


with  the  restrictions  on  variables  as  noted  in  the  assump- 
tions . 


Ill . Equipment  and  Procedures 


As  a secondary  purpose  of  this  study,  the  existing 
Mossbauer  spectroscopy  equipment  was  modified  to  simplify 
use  and  adjustment  and  to  improve  resolution.  This  chapter 
will  describe  the  equipment  used  in  data  acquisition,  as  well 
as  the  modifications  to  that  equipment.  The  final  sections 
of  this  chapter  will  describe  the  experimental  methods  used 
and  the  conditions  of  the  data  runs . 


The  major  components  of  the  equipment  used  were  manu- 
factured by  Ranger  Electronics,  and  consist  of  a velocity 
transducer,  a linear  amplifier/single  channel  analyzer,  a 
krypton-filled  proportional  counter,  and  a Mossbauer  control 
unit.  Additional  components  consisted  of  a RIDL  400  channel 
analyzer,  a Wavetek  oscillator,  a Ranger  temperature  control- 
ler, a heater,  a gold  (0.0?^  Fe)  - chromel  thermocouple,  and 
a digital  volt  meter.  These  components  are  those  used  by 
Skluzacek.  In  his  work  in  1976,  he  gives  an  overall  descrip- 
tion of  the  equipment  which  will  not  be  repeated  here  (Ref  10; 
5-11) . 


Cryostat  and  Drive  Mechanism.  A cryostat,  manufactured 
by  Janis,  Inc,  was  used  throughout  the  study  with  the  sample 
absorber  held  fixed  in  a combination  heater  and  mount . The 
heater/mount  and  moving  source  were  located  near  the  bottom 
of  the  cryostat,  and  were  attached  to  the  Mossbauer  drive 
unit  by  a 29  inch  support  and  an  extension  shaft.  The 
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possible  problem  of  a decrease  in  the  rate  of  change  of  vel- 
ocity near  zero  velocity,  as  reported  by  Skluzacek,  and  a need 
to  simplify  the  mounting  and  adjustment  of  the  drive  unit,  led 
to  the  following  modifications. 

The  heater/mount  support  and  motor  mount  were  changed  to 
the  configuration  indicated  in  Fig.  1.  The  changes  provided 
access  to  the  rods  forming  the  support  for  adjustment,  with- 
out removal  of  the  motor.  These  modifications  had  been  init- 
iated prior  to  the  start  of  this  study  by  Dr.  George  John, 
and  were  completed  by  the  writer. 

During  reassembly  and  testing  of  the  modified  support 
and  motor  mount,  two  problems  were  noted:  (1)  adjustment  of 
the  motor  extension  shaft  still  required  removal  of  the  motor 
and/or  the  Mossbauer  source;  and  (2)  operation  of  the  motor 
resulted  in  a lateral  vibration  of  the  extension  shaft.  It 
was  felt  that  the  vibration  was  due  to  the  flexibility  of  the 
small  diameter  tube  used  as  the  extension  shaft.  The  exten- 
sion shaft  shown  in  Fig.  2 was  designed  to  replace  the  exist- 
ing rod. 

The  new  shaft  was  constructed  of  thin  wall  l/4  inch  dia- 
meter stainless  steel  tubing.  The  adjustment  sleeve  was  made 
of  aluminum  to  reduce  weight,  and  the  threaded  inserts  were 
formed  from  steel  bolt  stock  center-bored  to  reduce  weight. 

The  source  mount  was  formed  from  brass  rod.  The  new  shaft 
displayed  (visually)  none  of  the  tendency  to  vibrate  dis- 
played by  the  original. • 

Moire  Gratinp;.  Skluzacek  also  reported  significant 
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Fig.  1 Modified  Motor  Mount  and  Support 


13. 


problems  with  non-linear  velocity  in  the  near  zero  region, 
due  to  the  use  of  a high  speed  Moire  grating  in  the  velocity 
measuring  system  (Ref  10).  Prior  to  the  start  of  this  study, 
the  grating  was  replaced  with  one  suitable  for  the  velocity 
range  to  be  used  (a  maximum  of  + 10  mm/sec) . The  velocity 
fits  obtained  with  the  new  grating  and  modified  extension 
shaft  were  typically  of  + 0.02  mm/sec,  and  were  obtained  with 
a first  order  polynomial  fit.  A calibration  run  made  with 
natural  iron  is  described  in  Chapter  IV. 

Temperature  Control.  Data  runs  were  accomplished  at 
both  room  and  liquid  nitrogen  temperatures.  Significant  tem- 
perature fluctuations  were  noted  at  both  temperatures . Room 
temperature  runs  were  accomplished  in  the  cryostat  with  the 
vacuum  system  operating,  but  large  fluctuations  in  the  lab- 
oratory ambient  temperature  resulted  in  a slow  drift  of  sample 
temperature  through  a range  of  + 5°C . 

Runs  at  liquid  nitrogen  temperature  were  initially  at- 
tempted at  a temperature  five  degrees  above  the  liquid  nitro- 
gen reference  used  for  the  gold-chromel  thermocouple.  This 
was  done  to  avoid  fluctuations  in  temperature  while  the  liquid 
nitrogen  was  being  serviced,  which  had  been  reported  by 
Skluzacek  (Ref  10:9)-  During  this  study,  runs  with  temper- 
ature fluctuations  of  up  to  three  degrees  Kelvin  were  noted; 
however,  the  subsequent  failure  of  the  temperature  indicating 
system  suggests  that  the  fluctuations  may  have  been  the  result 
of  the  malfunctioning  of  that  system,  and  not  actual  variations 
in  the  temperature  of  the  sample . 
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Source  and  Absorbers 


The  source  used  throughout  this  study  was  cobalt-57  in 
rhodium,  prepared  by  Amersham-Searle . The  source  was  origin- 
ally of  50  millicurie  strength,  but  had  decayed  to  approx- 
imately 16  millicuries  at  the  time  of  this  study. 

The  absorbers  were  manufactured  by  the  Batelle  Labor- 
atories at  Columbus,  Ohio,  and  were  made  available  through 
Dr.  Harold  Gegel,  of  the  Air  Force  Materials  Laboratory.  A 
list  of  the  samples  examined  is  given  in  Chapter  V . Each  of 
the  samples  was  manufactured  by  spin  cooling  a thin  stream  of 
the  melt  on  a cooled  rotating  drum.  All  of  the  samples  were 
similar  in  form,  consisting  of  thin,  narrow  strips.  Typically 
these  strips  were  0.001  inches  in  thickness  and  0.75  nim  in 
widths  The  thickness  showed  little  variation  from  sample  to 
samples  however,  the  width  of  a given  sample  would  range  from 
0.50  mm  to  1.2  mm.  Visual  examination  of  the  samples  wi '.h  a 
binocular  microscope  showed  the  edges  to  be  very  ragged  and 
the  surface  to  have  an  irregular  pattern  generally  running 
the  length  of  the  sample . 

A simple  calculation  was  made  to  determine  how  well  the 
samples  matched  the  conditions  developed  by  Shimony  for  max- 
imum single  line  absorption  (Ref  11).  While  the  samples  were 
not  expected  to  be  single  line  absorbers,  and  an  optimum  thick- 
ness would  exist  for  each  of  the  absorption  peaks,  the  cal- 
culation of  the  single  line  thickness  would  give  an  approx- 
imation of  the  actual  thickness  desired.  From  Shimony 's 
methods  and  tables,  the  calculated  thickness  was  0.001  inches. 
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Based  on  this  result,  the  samples  used  in  this  study  were 
formed  by  using  a single  thickness  of  the  material. 

Twenty  3/4  inch  strips  of  the  material  to  be  analyzed 
were  placed  side  by  side  on  ordinary  transparent  tape,  with 
the  strips  being  placed  on  the  tape  in  the  same  orientation 
that  they  were  cut  from  the  original  strip.  The  samples  were 
formed  large  enough  that  the  cut,  and  possible  distorted,  ends 
would  not  be  in  the  zone  used  for  analysis.  Similarly,  any 
kinked  or  distorted  portions  of  the  original  strip  were  not 
used. 

In  addition  to  the  amorphous  alloys,  two  other  absorbers 
were  used  for  calibration  of  the  system  and  testing  of  the 
curve  fitting  program:  a National  Bureau  of  Standards  sodium 
nitroprusside  crystal,  and  a natural  iron  foil. 

Data  Collection 

As  previously  noted,  each  sample  was  examined  at  botn 
room  and  liquid  nitrogen  temperatures.  For  each  of  these 
runs,  the  sample  was  mounted  between  two  0.75  ^ thick  disks 
of  boron  nitride.  This  assembly  was  mounted  approximately 
1.8  cm  from  the  midpoint  of  the  moving  source.  After  pass- 
ing through  the  absorber  and  disks,  the  radiation  passed 
through  two  mylar  windows  in  the  cryostat,  and  entered  the 
detector  through  a beryllium  window.  The  total  absorber  to 
detector  distance  was  18  cm.  The  arrangement  is  shown  in 
Fig.  3. 

Collection  times  for  the  various  runs  ranged  from  25  to 
6?  hours,  with  an  average  collection  time  of  50  hours.  Ex- 
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Fig.  3 Experimental  Arrangement 


perimental  conditions  during  the  runs  were  not  consistent.  A 
malfunctioning  heating  system  in  the  laboratory  caused  am- 
bient temperatures  to  vary  as  much  as  150  c during  the  runs. 

Several  sources  of  vibration  were  present  during  the 
runs.  The  vacuum  system  used  consisted  of  a fore  pump  and 
an  oil  diffusion  pump,  both  of  which  were  operated  contin- 
uously during  the  data  runs . The  vibration  generated  by 
these  pumps  was  transmitted  to  the  cryostat  through  the  re- 
quired connections.  Both  of  the  connections  used  v/ere  rig-id- 
ly  fixed  prior  to  reaching  the  cryostat,  but  vibration  could 
still  be  felt  when  the  pumps  were  operating. 

It  can  be  assumed  that  both  the  temperature  fluctuations 
and  the  vibrations  contributed  to  loss  of  intensity  and  to 
line  broadening.  The  extent  of  these  effects  has  not  been 
evaluated.  It  can  be  assumed  that  the  contribution  of  vibra- 
tion will  be  relatively  consistent  throughout  the  runs;  how- 
ever, the  effects  of  temperature  changes  may  be  both  material 
and  site  dependent,  and  are  therefore  much  more  difficult  to 
account  for  in  the  spectra. 

Data  Processing 

Raw  data  collected  and  stored  in  the  multi-channel  ana- 
lyzer consisted  of  time-base  oscillator  counts  in  channels  0 
to  15.  velocity  data  in  channel  8 and  every  tenth  channel 
thereafter,  and  the  obtained  spectrum  in  the  remaining  chan- 
nels. After  completion  of  a data  run,  the  counts  in  the  first 
70  channels  were  recorded  by  an  automatic  typewriter  for  run 
identification.  The  entire  spectrum  was  then  recorded  on 
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paper  tape  with  a Tally  paper  tape  punch.  The  paper  tape  was 
converted  to  Hollerith  punched  computer  cards  by  use  of  a 
permanent  library  program.  Initial  processing  consisted  of 
the  calculation  of  the  velocity  curve  and  the  generation  of 
a counts  vs.  velocity  plot.  Further  processing  of  the  data 
was  made  using  the  GENFIT  program,  and  is  discussed  in 
Chapter  IV . 


20 


IV . Data  Analysis 


The  experimentally  obtained  data  was  analyzed  using  a 
least-squares  minimization  curve-fitting  program  originally 
developed  by  the  Argonne  National  Laboratory.  This  program, 
named  GENFIT,  requires  a user-supplied  subroutine  to  describe 
the  spectrum  being  analyzed.  The  GENFIT  program  varies  those 
parameters  specified  by  the  user  to  obtain  the  best  least- 
squares  fit.  The  parameters  used  in  this  study  are  those  de- 
scribed in  Chapter  II,  which  determine  the  position,  width, 
and  height  of  the  absorption  peaks  generated  by  each  combin- 
ation of  atoms . Numerous  changes  to  the  GENFIT  program  and 
the  READIN  subroutine  were  made  to  simplify  use  and  to  elim- 
inate those  portions  not  applicable  to  Mossbauer  spectroscopy. 
These  changes  are  described  in  Appendix  A . 

The  user-supplied  subroutine,  called  CALFUN,  provides 
the  means  by  which  the  analyst  implements  an  assumed  model. 
This  chapter  will  discuss  the  program  GENFIT,  the  subroutine 
CALFUN,  the  criteria  used  for  goodness-of-f it , and  finally, 
the  techniques  used  in  the  employment  of  the  program. 

Program  GENFIT 

The  program  GENFIT  was  modified  during  this  study  to  add 
a line-printer  plotting  capability  and,  as  noted,  to  simplify 
its  use.  The  portions  of  the  program  performing  the  actual 
minimization  were  left  essentially  unchanged.  The  subroutine 
READIN,  written  by  Skluzacek,  was  modified  to  accommodate  a 
change  in  the  channel  numbers  used  for  multiplexed  velocity 
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data,  and  to  allow  for  partial  overflows  in  the  data  spectrum. 

A listing  of  the  GENFIT  main  program  and  the  subroutine 
I READIN  is  given  in  Appendix  A . 

Subroutine  CALFUN 

The  subroutine  CALFUN  is  called  by  the  minimization 
routine  to  describe  the  spectrum  being  calculated.  To  reduce 
the  execution  time  of  the  program,  the  simplified  Hamiltonian, 
as  well  as  the  binomial  probabilities,  were  reduced  to  numer- 
ical values,  and  used  as  constants  ’-vherever  possible. 

The  subroutine  used  a total  of  1?  variables  to  Implement 
the  model  described  in  Chapter  II.  The  required  variables  i! 

were: 

1.  Baseline:  the  average  number  of  counts  in  the  back- 
ground spectrum. 

;•  2.  Magnetic  field:  six  va‘'.ues  describing  the  intern.?! 

magnetic  field  for  each  of  the  combinations. 

3.  Total  intensity:  total  absorption  intensity  divided  : 

1 

by  TT  . i 

' 4.  Isomer  shift:  self  explanatory. 

5.  Linewidth:  three  values  describing  the  full  width  at 
half  maximum,  one  each  for  peaks  1-6,  2-5i  and  3“^* 

6.  Ratios:  five  values  representing  the  ratio  of  a 
given  peak  to  the  height  of  peak  six. 

A listing  of  the  subroutine  CALFUN  is  included  in  Appendix  B. 

Goodness-of-Fit 

Two  criteria  were  used  to  measure  the  goodness-of-f it 
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obtained  by  a given  computer  run.  One  of  these  was  objective 
(Chi-squared),  and  the  other  subjective  (the  visual  appear- 
ance of  the  calculated  spectrum) . These  criteria  were  pri- 
marily used  to  compare  the  results  of  various  computer  runs, 
and  only  provide  a secondary  indication  of  the  validity  of 
the  model. 

A value  of  Chi-squared  was  generated  by  the  GENFIT  pro- 
gram and  differs  slightly  in  definition  from  that  normally 
used.  Because  the  total  number  of  data  points  considered  in 
each  of  the  spectra  varies  slightly,  the  value  of  Chi-squared 
is  defined  by 


Chi2= 


(data  pointi-calculated  pointi)2 
data  pointi 


(5) 


The  subjective  criterion  of  visual  appearance  was  re- 
quired because  of  the  program  GENFIT 's  tendency  to  "lock  in" 
to  local  minima,  rather  than  to  absolute  minima.  In  these 
cases,  the  calculated  spectrum  could  have  a reasonable  value 
of  Chi-squared,  and  yet  bear  little  resemblance  to  the  exper- 
imental data.  In  such  cases,  the  visual  appearance  of  the 
generated  plots  was  the  controlling  criterion. 


Use  of  the  Program 

The  use  of  a model  which  forced  a binomial  distribution 
on  the  contributions  of  each  combination  allowed  a simpler 
means  of  analysis  than  would  have  been  possible  with  free 
variation  of  the  contributions.  The  plot  of  counts  vs.  vel- 
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ocity  obtained  in  initial  processing  was  used  to  determine 
approximate  values  for  the  various  parameters.  For  example, 
the  site  or  combination  having  the  greatest  contribution 
would  be  assumed  to  have  an  internal  magnetic  field  such  that 
its  peaks  one  and  six  would  fall  directly  under  peaks  one  and 
six  of  the  experimental  peaks.  Using  similar  methods,  and 
assuming  a distribution  of  values  proportional  to  the  percent 
of  iron  present  at  a given  site,  a reasonable  estimate  for 
each  of  the  parameters  can  be  found. 

The  initial  run  of  the  program  GENFIT  was  made  using  a 
very  limited  number  of  iterations,  and  without  requiring  con- 
vergence. The  results  of  this  run  allow  the  detection  of  any 
gross  errors  in  the  initial  estimates,  as  well  as  disclosing 
any  singularities  in  the  values.  It  was  noted  in  the  data 
analysis  that  certain  combinations  of  values  "trap"  the  min- 
imization scheme  in  a non-converging  iteration  of  values . 

Once  the  initial  values  have  been  confirmed  or  modified 
by  these  short  runs,  sucessive  runs  can  be  made  using  the 
results  of  the  previous  run.  Two  exceptions  exist  to  this 
general  rule.  The  first  is  the  use  of  the  variable  called 
baseline.  The  minimization  scheme  treats  baseline  as  a free 
variable,  although  it  represents  the  value  of  the  average 
number  of  counts  in  the  background  spectra.  It  is  recom- 
mended that  this  value  be  calculated  manually,  and  that  this 
value  be  used  for  all  computer  runs  of  that  spectrum.  The  sec- 
ond exception  occurs  when  the  calculated  spectrum  deviates 
significantly  from  the  experimental.  Normally,  in  these 
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cases,  the  program  has  "locked  in"  to  a local  minimum,  and  it 
is  necessary  to  adjust  manually  one  or  more  of  the  parameters 
away  from  the  program  generated  value  and  re-execute  the  pro- 
gram. 

Further  comments  on  the  use  of  the  program  GENFIT  can  be 


found  in  Ref  10. 


V . Results 


The  Mossbauer  pattern  obtained  for  each  of  the  materials 
analyzed  and  the  best-fit  calculated  spectrum  are  shown  in 
Figs.  4 to  17  in  this  chapter.  Each  of  the  calculated  spectra 
were  generated  using  the  model  described  in  Chapters  II  and 
IV.  Table  I presents  the  values  obtained  for  the  1?  para- 
meters used  in  the  curve  fitting,  as  well  as  the  conditions 
of  each  run. 

Runs  are  identified  by  sample  composition,  followed  by 
a run  number.  For  example,  FE80B20  #2  would  be  the  second 
run  with  a sample  composed  of  80??  iron  and  20%  boron. 

Certainty 

Uncertainty  values  for  the  parameters  listed  in  Table  I 
are  not  listed,  in  order  to  reduce  confusion  in  the  table.  In 
addition,  the  magnitude  of  the  i.ncertainties  is  primarily  a 
function  of  the  computer  time  allowed  for  the  program  to  con- 
verge to  a minimum  least  squares.  A list  of  the  typical 
values  of  uncertainty  for  each  of  the  variables  is  included 
here,  and  represents  the  worst  values  .found. 

For  internal  magnetic  fields,  the  value  was  typically 
+0.5  kOe  for  the  five  largest  fields;  for  the  smallest  cal- 
culated field,  the  value  was  typically  +2.5  kOe.  For  line- 
widths,  the  value  was  typically  0.01  mm/sec  or  less.  Isomer 
shift  values  were  + 0.001  mm/sec,  and  values  for  the  peak 
height  ratios  were  + 0.02  or  less.  Other  factors  relating 
to  the  certainty  of  the  calculated  values  are  discussed  in 
Chapter  VI . 


Table  I 


Table  One  summarizes  the  results  of  the  calculated 
spectra.  Abbreviations  used  in  the  table  are  Temp  for 
temperature  in  degrees  Kelvin,  Hfld  for  internal  magnetic 
field  expressed  in  kOe,  Chi  Sqd  for  Chi-squared,  LW  for  line- 
width,  which  is  followed  by  the  applicable  peak  numbers.  Iso 
Sft  for  isomer  shift,  and  Rat  for  the  ratio  of  the  indicated 
peak  height  to  peak  number  six.  Throughout  this  report, 
peaks  are  numbered  from  the  peak  with  the  most  negative 
velocity  as  number  one,  to  the  peak  with  the  most  positive 
velocity  as  peak  number  six. 
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Measured  and  Calculated  Parameter  Values 
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Symbols  and  units  defined  on  page  27- 


Table  I (Continued) 


Linewidths  in  mm/sec. 
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Fig.  7 Mossbauer  Spectrum  of  FeonBon  at  75°K 
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•Fig.  11  Mossbauer  Spectrum  of  Fe3QP3^^Cj[j,^^By  at  75°K 
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VI.  Conclusions  and  Recommendations 

The  results  obtained  in  this  study  lead  to  several  con- 
clusions. While  the  results  indicate  a general  failure  of 
the  model  in  describing  the  structure  of  the  samples,  two 
notable  exceptions  will  be  discussed  in  this  chapter.  In 
addition,  several  conclusions  independent  of  the  model  can 
be  drawn.  This  chapter  will  first  discuss  the  shortcomings 
of  the  model  used,  the  two  materials  which  appear  to  be  ex- 
ceptions, the  conclusions  independent  of  the  model,  and, 
finally,  will  offer  several  recommendations  toward  further 
study  of  these  amorphous  materials. 

Structural  Model 

The  use  of  the  substitutional-solid  model  produced  cal- 
culated spectra  that  were  generally  satisfactory  visually. 

However,  the  parameter  values  generated  in  the  curve  fitting 
raise  a question  as  to  the  validity  of  the  results. 

The  values  of  Chi-squared  obtained  ranged  from  443  for 
the  "best”  fit  to  3»^99  for  the  "worst".  If  the  model  were 
to  describe  all  of  the  materials  adequately,  the  values  of 
Chi-squared  obtained  would  be  expected  to  have  a random  dis- 
tribution over  a much  smaller  range,  in  the  limit  approaching 
the  number  of  data  points  fitted.  The  room  temperature 
spectra  of  the  parent  material  approached  this  value  with  a 
Chi-squared  of  443  and  approximately  3^0  data  points;  how- 
ever, the  same  material  at  yielded  a Chi-squared  value 

of  1,187.  This  trend  was  repeated  for  all  except  one  of  the 
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spectra  obtained.  The  generally  high  values  of  Chi-squared 
and  the  non-random  distribution  of  values  with  changed  tem- 
perature indicate  a failure  of  the  model  to  describe  ade- 
quately the  structure  of  all  the  samples . The  assumption 
that  constituent  atoms  could  be  treated  as  iron  and  non- 
iron is,  as  anticipated,  invalid. 

Significant  here,  however,  are  the  results  obtained  with 
FegQB2o*  Here  the  constituents  are  iron  and  non- iron,  and 
the  probability  of  a site  being  occupied  by  either  iron  or 
boron  should  follow  a binomial  distribution.  Yet  the  results 
obtained  were  not  the  best  obtained;  in  fact,  these  were  the 
highest  values  of  Chi-squared  obtained.  This  discrepancy 
could  be  explained  if  the  iron-boron  compound  acted  in  the 
manner  described  by  Fujita  for  a body  centered  cubic  iron- 
rich  compound  (for  which  second  neighbor  effects  are  not 
negligible) , while  the  more  comjlex  compounds  formed  not 
only  the  body  centered  cubic,  but  also  contained  a high  per- 
centage of  randomly  distributed  interstitial  atoms.  The  com- 
bined effects  of  an  absorber's  first  neighbors  and  closest 
interstitial  neighbors  might  be  more  nearly  described  by 
this  study's  model  than  is  the  binary  compound  itself. 

The  values  of  linewidth  and  peak  ratio  generated  in  the 
curve  fitting  process  results  in  a second  problem.  The 
values  obtained  for  linewidths  are  similar  to  those  obtained 
by  Tsuei,  ranging  from  I.03  mm/sec  to  0.37  mm/sec  (Ref  8:605)* 

Fujita,  using  his  model,  reports  values  for  linewidths  "a 

few  tens  of  percent  broader"  than  that  obtained  for  pure  iron 

I 
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(Ref  6:232).  In  a calibration  run  made  with  natural  iron, 
an  average  linewidth  of  0.24  mm/sec  was  obtained,  indicating 
that  the  experimentally  obtained  linewidths  are  considerably 
in  excess  of  a few  tens  of  percent  larger  than  the  natural 
linewidths.  The  implication  is  that  the  values  of  linewidth 
were  increased  by  the  curve  fitting  routine  in  order  to 
"fill"  the  absorption  peaks,  and  that  contribution  from  more 
site  combinations  was  required.  It  is  not  felt  that  using 
more  than  six  of  the  twelve  sites  in  this  model  would  have 
improved  the  values  of  linewidth.  More  combinations  with 
nearly  equal  contributions  would  be  required. 

The  relative  peak  heights  obtained  also  indicate  a pos- 
sible problem.  If  the  spins  are  randomly  oriented  in  the 
sample,  the  value  obtained  for  peak  height  should  be  in  an 
approximate  ratio  of  3:2: 1:1:2; 3*  This  set  of  ratios  was  not 
obtained  with  the  experimental  data.  The  obtained  ratios 
can  indicate  a non-random  distribution  of  spins.  However, 
again  the  question  of  whether  enough  combinations  were  con- 
sidered to  properly  describe  the  spectrum  must  be  considered. 
The  ratios  obtained  may  not  be  significant. 

Evidence  exists  to  suggest  that  the  results  obtained 
are  possibly  a function  of  both  of  these  causes.  The  ex- 
perimental spectra  were  not  perfectly  symmetrical,  and  the 
ratios  changed  significantly  with  temperature:  both  suggest- 
ing an  other-than-random  spin  in  the  absorber  atoms.  On  the 
other  hand,  allowing  the  ratios  to  vary  arbitrarily,  as  well 
as  requiring  the  ratios  of  each  site  contribution  to  be  the 
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same,  places  a great  uncertainty  on  the  values  obtained.  The 
question  remains  unresolved. 

Exceptions 

In  those  samples  containing  boron,  the  percentage  of 
boron  increases,  and  the  amplitude  of  the  outermost  peaks 
shows  a nearly  proportional  increase  until  seven  percent 
boron  is  reached.  Above  seven  percent,  the  amplitude  de- 
creases. In  the  seven  percent  boron  sample,  the  outer  peaks 
have  an  intensity  O.OI5  greater  than  peaks  two  and  five.  The 
two  percent  silicon  sample  has  a similar  distribution  of  peak 
heights  with  the  outer  peaks  only  0.005  less  intense  than 
peaks  two  and  five.  It  is  only  in  these  two  samples  that 
peaks  2,  3i  and  5 show  a 1:1:1 :1  ratio.  These  effects  are 
only  evident  in  the  runs  at  liquid  nitrogen  temperature  (Figs. 
10,  14). 

Two  possible  causes  exist  to  explain  these  significant 
changes  in  the  samples  as  they  are  cooled.  The  first  is  that 
there  is  a significantly  different  vibrational  mode  of  the 
resonating  atoms , resulting  in  a change  in  the  relative  line 
intensities  of  the  contributing  sites.  A second  possibility 
is  a different  angular  dependence  of  tne  nuclear  Zeeman 
effect.  The  first  of  these  possibilities  suggests  a non- 
cubic  symmetry,  and  the  second  suggests  a material  with  a 
preferred  orientation  (Ref  12:16).  In  either  case,  the  re- 
sults indicate  a different  structure  for  these  two  materials 
at  liquid  nitrogen  temperature  than  at  room  temperature,  and 
as  compared  to  the  remaining  materials . 
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Model  Independent  Results 


The  patterns  obtained  through  Mossbauer  spectroscopy 
are  only  a measure  of  the  effect  that  surrounding  atoms  have 
on  the  absorber  nucleus,  and  it  is  conceivable  that  two  dif- 
ferent combinations  of  neighbor  atoms  could  result  in  the 
same  Mossbauer  spectrum.  Fujita  reports  obtaining  the  same 
spectra  for  three  different  (atomic  percentages)  iron-phos- 
phorous-carbon amorphous  alloys  (Ref  6:231) • In  this  study, 
the  room  temperature  spectra  of  FegoPn  7*^6. 3^ ^2  differ  sig- 
nificantly only  in  peak  heights.  As  previously  discussed,  the 
peak  height  ratios  generated  by  the  model  are  of  questionable 
value.  Similar  compositions  could  lead  to  similar  values  for 
the  internal  magnetic  fields,  and  to  different  preferential 
directions  of  spin,  or  angular  dependence  of  the  Zeeman  effect. 
The  different  values  of  isomer  shift  obtained  indicate  that 
the  two  spectra  are  unique,  but  further  study  is  i .dicated  to 
determine  if  these  two  materials  do  have  unique  spectra. 

While  the  model  has  failed  to  describe  adequately  the 
materials,  the  consistent  increase  in  the  value  of  Chi- 
squared  obtained  with  decreasing  temperature,  combined  with 
the  visual  change  in  the  spectra,  indicates  a change  in  the 
structure  of  the  materials  with  decreasing  temperature. 

The  wide  range  of  values  of  Chi-squared  also  implies 
differences  in  the  structures  of  the  various  materials,  even 
at  room  temperature.  A single  model  may  not  be  sufficient 
to  describe  all  of  the  materials. 


In  summary,  the  solid  solution  twelve  site  model  was 
inadequate  to  describe  fully  the  materials  analyzed.  This 
conclusion  is  primarily  based  on  the  values  of  Chi-squared 
and  linewidth  obtained  in  the  curve  fits.  Further  research 
is  indicated  to  determine  which  portions  of  the  model  need  to 
be  revised,  and  whether  the  two  exceptions  noted  can  be 
fully  described  by  the  model.  The  model  has  established 
significant  trends  in  the  various  spectra,  and  may  therefore 
be  a valuable  tool  in  establishing  the  initial  values  that 
would  be  required  in  using  a more  complex  model  of  the  amor- 
phous structure  of  these  materials.  Finally,  the  results  of 
this  study  indicate  that  the  model  used  may  be  temperature 
dependent  at  lower  temperatures,  as  Tsuei's  results  indicated 
at  above  room  temperature  (Ref  8) . 

Recommendations 

An  evaluation  of  these  materials  using  the  iron-rich 
body  centered  cubic  model  of  Fujita  would  be  the  next  log- 
ical step  in  the  analysis  of  these  amorphous  materials.  A 
jump  to  the  use  of  all  14  sites  is,  however,  not  recommended. 
A progression  in  the  number  of  sites  (and  interstitial  sites) 
considered  should  be  made  until  acceptable  values  of  line- 
width  are  obtained.  This  recommendation  is  based  on  a con- 
viction that  the  minimum  number  of  sites  considered  with 
consistent  results  will  be  the  most  reasonable  approximation 
to  the  true  structure  of  the  materials. 

The  second  recommendation  is  that  regardless  of  the 
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model  assume  in  future  studies,  spectra  should  be  obtained 
throughout  the  range  of  temperatures  from  liquid  nitrogen 
to  room  temperature.  It  remains  to  be  determined  whether 
the  changes  with  temperature  noted  in  this  study  are  gradual 
over  the  range,  or  represent  phase  changes  occurring  at  a 
specific  temperature. 

The  results  obtained  with  seven  percent  boron  and  two 
percent  silicon  require  confirmation.  The  remaining  samples 
of  parent  plus  silicon  (not  examined  in  this  study)  need  to 
be  evaluated  to  determine  if  they  exhibit  the  same  character- 
istics. In  addition,  confirmation  that  the  changes  noted  are 
material,  and  not  manufacture  dependent,  is  required.  Only 
one  sample  of  each  material  was  available  in  this  study. 

A final  recommendation  is  that  the  results  of  this  study 
be  compared  with  the  x-ray  diffraction  and  other  tests  now 
being  accomplished  by  the  Air  Force  Materials  Laboratory. 
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Appendix  A 

Program  GENFIT  and  Subroutine  READIN 
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Program  GENFIT  is  the  controlling  routine  for  the  least 
squares  fitting  program  used  in  this  study.  As  obtained, 
the  program  was  designed  for  either  Mossbauer  or  general 
curve  fitting  use.  To  reduce  execution  time  and  simplify 
use,  this  portion  of  the  program  was  significantly  rewritten. 
Similarly,  subroutine  READIN  was  modified  to  be  compatible 
with  the  revised  main  program,  and  was  also  modified  to 
accept  data  in  a new  format , 

The  input  for  both  of  these  routines  is  read  from  cards 
in  the  following  order  and  consists  of: 

1.  Two  cards  containing  the  spectrum  title  or  desired 
comments.  Each  of  these  cards  may  contain  up  to  80 
characters.  The  information  is  stored  in  array  TITL. 

2.  One  card  with  six  integer  values  separated  by 
commas  or  blank  spaces.  These  variables,  and  the 
order  in  which  they  should  appear,  are: 

a.  NEAR:  the  number  of  variables  to  be  used  in 
subroutine  CALFUN  to  describe  the  spectrum. 

b.  VMXl:  the  maximum  channel  number  of  the  left 
branch  of  the  velocity  spectra  prior  to  zero 
area  irregularities. 

c.  VMIN2:  minimum  channel  number  of  velocity 
spectra  in  right  branch  after  zero  area 
irregularities . 

d.  TOVR:  the  number  of  overflows  in  the  time  base 
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oscillator  counts. 

e.  IPLOT:  Control  parameter  for  the  type  of 
output  plot  desired: 

IPL0T=0  CALCOMP  plot  of  experimental  and 
calculated  spectra 

IPL0T=10  Line  printer  plot  of  experimental 
and  calculated  spectra 

IPLOT=any  other  value  CALCOMP  plot  of 
experimental  data,  calculated 
spectrum,  F-contribution , and 
residuals 

f.  MAXFUN:  The  maximum  number  of  calls  of  sub- 
routine CALFUN  desired. 

3.  NPAR  cards  with  the  name  of  the  variable  in  the  first 
six  columns,  and  the  initial  estimate  of  the  value 

in  columns  six  through  twelve.  These  cards  are  read 
in  an  A6e6  format  and  may  be  combined  onto  one  or 
two  cards  in  consecutive  fields  of  twelve  columns. 

4.  Forty  cards  containing  the  raw  data  as  generated  by 
established  tape  processing  program. 

5.  One  card  with  a variable  number  of  integer  values, 
read  in  free  field  format,  in  the  following  order: 

a.  ORDER:  The  order  of  the  polynomial  velocity 
fit  desired. 

b.  DOVR:  The  number  of  overflows  in  the  Mossbauer 
data.  If  the  data  has  not  completely  overflowed, 
the  higher  value  is  used. 
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c.  NZONES:  The  number  of  zones  to  be  used  in  the 
background  polynomial  fit. 

d.  ZONEL  and  ZONER:  The  start  and  stop  channel 
numbers  of  the  number  of  zones  specified  in  c. 

A listing  of  the  revised  GENFIT  main  program  and  the 
READIN  subroutine  are  included  in  this  appendix. 
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Appendix  B 


Subroutine  CALF UN 

Subroutine  CALFUN  is  used  to  describe  the  calculated 
spectrum  using  1?  variables,  adjusted  by  the  GENFIT  least- 
squares  minimization  routine.  This  version  of  CALFUN  uses 
the  simplified  Hamiltonian,  further  reduced  to  a numerical 
form  for  speed  of  execution.  In  order  to  calculate  the  peak 
positions,  the  routine  used  variables  of  magnetic  field  (one 
for  each  assumed  site),  and  isomer  shift.  To  compute  the 
Lorenzian  line  shapes,  the  routine  uses  the  binomial  prob- 
abilities multiplied  by  an  artificial  variable,  called  total 
intensity,  to  find  individual  site  and  peak  intensities.  To 
implement  the  model,  linewidths  are  variable,  with  the  re- 
striction that  the  nth  peak  of  each  site  will  have  the  same 
linewidth.  The  factors  called  ratio  [[RAT(x)]  allow  the 
individual  peak  heights  for  a given  site  to  assume  whatever 
ratio  yields  the  best  fit — again  with  the  restriction  that  the 
nth  peak  of  each  site  be  in  the  same  ratio.  Peak  heights  are 
normalized  to  peak  six  (the  peak  with  the  greatest  positive 
velocity) . 

The  final  spectrum  is  the  sum  of  the  six  sites'  contrib- 
utions, each  consisting  of  six  Lorenzian  profiles  centered  at 
the  velocities  determined  by  the  Hamiltonian. 

A listing  of  the  subroutine  CALFUN  is  included  in  this 
appendix. 
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